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Elliptic Flow of Protons and Antiprotons in Au+Au Collisions 
at -^snn = 7.7-62.4 GeV within Alternative Scenarios of Three-Fluid Dynamics 

Yu.B. Ivanov^'B 

^ Kurchatov Institute, Moscow RU-123182, Russia 

Analysis of elliptic flow of protons and antiprotons in Au+Au collisions is performed in a wide 
range of incident energies y/sNN = 7.7-62.4 GeV. Simulations has been done within the three-fluid 
model employing a purely hadronic equation of state (EoS) and two versions of the EoS involving 
deconfinement transition: an EoS with the first-order phase transition and that with a smooth 
crossover transition. It is found that the proton data are reproduced approximately to the same 
extent within all of the scenarios, including the hadronic one, while the deconfinement scenarios 
look certainly preferable for the antiproton elliptic flow. The fact that difference between elliptic 
flows of protons and antiprotons decreases with the incident energy rise is a consequence of reducing 
baryon stopping rather than an onset of deconfinement. 
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I. INTRODUCTION 

The large azimuthal anisotropic flow at the Relativis- 
tic Heavy Ion ColHder (RHIC) is believed to be a con- 
clusive evidence for the creation of dense partonic mat- 
ter in ultra-relativistic nucleus- nucleus collisions. This 
anisotropy is described by flow parameters defined as the 
proper Fourier coefficients Vn of the particle distributions 
in azimuthal angle with respect to the reaction plane [1| . 
The major bulk of research both theoretical and experi- 
mental has been done on the elliptic flow ('^2). 

Recently new data of the STAR Collaboration [2] on 
transverse-momentum dependence of the elliptic flow of 
identified particles in the incident energy range of y^SAT/v 
= 7.7-62.4 GeV were reported. These data were taken 
within the Beam Energy Scan (BES) program proposed 
at RHIC. The BES program was proposed in order to 
study possible onset of deconfinement, as well as possibly 
to identify the location of the critical end point that ter- 
minates the crossover transition at small quark-chemical 
potential to a first-order phase transition at higher quark- 
chemical potential [3|]. 

This work aims to analyze the new STAR data [2l] 
within different scenarios (with and without deconfine- 
ment) of heavy- ion collisions in order to draw conclusions 
on the matter produced in these collisions. The simula- 
tions were performed within a model of the three-fluid 
dynamics (3FD) (^ employing three different equations 
of state (EoS): a purely hadronic EoS [5] (hadr. EoS) and 
two versions of EoS involving the deconfinement transi- 
tion [6|. These two versions are an EoS with the first- 
order phase transition (2-phase EoS) and that with a 
smooth crossover transition (crossover EoS). Details of 
these calculations are described in Ref. [3]. First results 
of the 3FD simulations within alternative scenarios of 
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heavy- ion collisions are reported in Refs. [3-[9| dedicated 
to analysis of the baryon stopping and particle produc- 
tion. Results on transverse flow in the energy range of 
AGS (Alternating Gradient Synchrotron) and SPS (Su- 
per Proton Synchrotron)but only within hadronic EoS 
were presented in Refs. pjOl, Hlj. 

The 3FD model [3] treats the process of a nuclear col- 
lision from the very beginning, i.e. from the stage of 
incident cold nuclei, to the final stage of freeze-out. Con- 
trary to the conventional hydrodynamics, where local in- 
stantaneous stopping of projectile and target matter is 
assumed, a specific feature of the 3FD is a finite stop- 
ping power resulting in a counter-streaming regime of 
leading baryon-rich matter. The basic idea of a 3-fluid 
approximation to heavy-ion collisions |12l-[l^ is that at 
each space-time point a generally nonequilibrium distri- 
bution of baryon-rich matter can be represented as a sum 
of two distinct contributions initially associated with con- 
stituent nucleons of the projectile (p) and target (t) nu- 
clei. In addition, newly produced particles, populating 
the mid-rapidity region, are associated with a fireball (f) 
fluid. Therefore, the 3-fluid approximation is a minimal 
way to simulate the finite stopping power at high incident 
energies. 

The main observation of Ref. [2| is a strong differ- 
ence in V2{pt) between particles and their corresponding 
antiparticles. It is argued that it cannot be explained 
in a purely hydrodynamic approach since particles and 
anti-particles have the same mass. The evolution of the 
elliptic flow in the transient BES energy range was ad- 
dressed in a number of theoretical works |15|-[2Q| . 

In Refs. [15|, [l6| the difference between particles and 
antiparticles is explained by the presence of a vector 
mean field potential which is repulsive for particles and 
attractive for antiparticles. In [l5| it is introduced at 
the quark level within the Nambu-Jona-Lasinio (NJL) 
model, while in 'l6y, at the hadronic level. As a conse- 
quence of these potentials, antiparticles are attracted by 
the matter and are trapped in the system, whereas parti- 
cles feel a repulsive force and have the tendency to leave 



the system along the participant plane. With the poten- 
tials included, a fair qualitative agreement was achieved 

In Ref. pJ], a hybrid (hydrodynamical plu s Ultra- 
relativistic Quantum Molecular Dynamics [21|) calcu- 
lation was performed. The effect for the protons- 
antiprotons primarily results from the annihilation pro- 
cess. Antiprotons moving in the out-of-plane direction 
encounter more protons to annihilate with than those 
moving in the in-plane direction. Another effect dis- 
cussed in this paper is related to the event plane calcula- 
tion. It was claimed that fluctuations in this calculation 
can bias the event plane to be rotated towards the most 
abundantly produced particles. This would, for exam- 
ple, increase the V2 values for protons and reduce them 
for anti-protons. 

The 3FD model proposes a more plausible (as com- 
pared to Refs. [15|, [l6|) explanation of the difference be- 
tween particle and antiparticle elliptic flow in terms of 
three interacting fluids: the fireball, consisting of parti- 
cles newly produced near the spacial center of the collid- 
ing system, and two baryon-rich fluid {p and t) associated 
with leading particles which traversed the whole system 
and are finally located in longitudinally peripheral re- 
gions. The particle-antiparticle difference explained by 
the increase of nuclear stopping in heavy ion collisions 
with decreasing energy. When the nuclear stopping be- 
comes strong, the mid-rapidity quantities are determined 
not only by particles newly produced near the spacial 
center (the f-fluid) but also contributed by leading par- 
ticles (the p- and t-fluids). The center and peripheral 
regions differently contribute to the mid-rapidity elliptic 
flow of different species, because they have different con- 
tent of particles-antiparticles (quarks-antiquarks). This 
naturally results in different V2 of particles and antipar- 
ticles because the center and peripheral regions have dif- 
ferent V2 patterns. This explanation is similar in some 
features to those proposed in Refs. [Tsl, [l9|. Contrary 
to Refs. [Isl, [IqI, the 3FD model calculates contribution 
from peripheral regions, rather than does assumptions on 
them, and does not employs the quark coalescence. 



II. COMPARISON WITH DATA 

The elliptic flow is proportional to the spatial 
anisotropy [22|, l23|. Usually, for this purpose one uses 
the eccentricity e defined by 



usual Woods-Saxon profile of nuclear density 



(y2) + (x2) 



(1) 



Mean values of spacial transverse coordinates (x^) 
(out of the reaction plane) and {y'^) (in the reaction 
plane) are usually calculated with either the wounded- 
nucleon (WN) or the binary-collision (BC) weights, for 
details see Ref. UM • These calculations are based on the 
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where po is the normal nuclear density, Ra = 1.12A^^^ 
is the raduis of a nucleus with mass number A, and d is 
a diffuseness of the nuclear surface. 

As long as the eccentricity is small, elliptic flow should 
be directly proportional to the eccentricity. For numeri- 
cally large eccentricities the direct proportionality could 
break in principle, but as was shown in the very first 
hydrodynamic calculation by Ollitrault [22| the propor- 
tionality holds well even for rather large values of e. 
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FIG. 1: Spatial eccentricity e as a function of impact 
parameter in Au+Au collisions for different surface diffuse- 
nesses of the Au nucleus (d) and different weights of aver- 
aging: the wounded-nucleon (WN) and the binary-collision 
(BC) weights [IJ. 

Within the 3FD model the initial nuclei are represented 
by sharp-edged spheres, i.e. the zero diffuseness (d = 0). 
This is done for stability of the incident nuclei before 
collision. This circumstance essentially affects the eccen- 
tricity. The results obtained with d = and the reahs- 
tic value of d = 0.6 fm calculated with BC weights are 
shown in Fig. [TJ As seen, the (d = 0)-result noticeably 
exceeds the eccentricity for the physically realistic value 
of d = 0.6 fm. Moreover, the (d = 0.6 fm)-result with BN 
weights practically coincides with the eccentricity calcu- 
lated with WN weights. The latter is considered as a 
realistic eccentricity and is accepted in the experimental 
analysis Ref. Q • 

The overestimation of e in the 3FD model naturally re- 
sults in the respective overestimation of the elliptic flow. 
There are two ways to compensate this overestimation: 
either the impact parameter of the collision should be re- 
duced to get a reasonable value of e for a considered cen- 
trality bin or the calculated value of V2 should be rescaled 
with the factor of eBN{d = 0.6 fin.)/eBN{d = 0), because 
V2 cx £, as mentioned above. The latter method is ap- 
plied in this paper, i.e. all the V2 values displayed below 
are rescaled with the above factor. As the calculations 
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FIG. 2: Elliptic flow of protons and antiprotons at mid-rapidity in central (0%-10%) Au+Au collisions at ^'snn = 7.7, 
11.5, 19.6, 27, 39 and 62.4 GeV as a function of transverse momentum. The 3FD calculations are performed at b = 3 fm. 
Experimental data are from Star Collaboration [2]. 
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FIG. 3: The same as in Fig. |2]but for mid-central (10%-40%) Au+Au collisions. The 3FD calculations are performed at b - 
6 fm. 
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FIG. 4: The same as in Fig. [2] but for peripheral (40%-80%) Au+Au colhsions. The 3FD calculations are performed at b - 
fm. 



are performed at fixed impact parameters, only proton 
and antiproton elliptic flow is considered for which data 
with comparatively narrow centrality selection are avail- 
able 0. 

Results of calculations of proton and antiproton ellip- 
tic flow at mid-rapidity and their comparison with STAR 
data [21 for three centrality bins are presented in Figs. [21 
[3] and [H The quality of data reproduction is quite rea- 
sonable for calculations at fixed impact parameters that 
are representative for the centrality bins considered. It 
should be kept in mind that in spite of the best repro- 
duction at the top considered energy of 62.4 GeV, results 
at this energy are still not quite accurate, since an ac- 
curate computation requires unreasonably high memory 
and CPU time. Precisely due to this reason the calcula- 
tions of peripheral (6 = 8 fm) collisions at 62.4 GeV have 
not been done. 

All scenarios, including the purely hadronic one, quan- 
titatively reproduce the proton data approximately to the 
same extent. It simply means that all scenarios represent 
strongly collective (fluid-like) behavior of the system that 
results in a large elliptic flow. The value of V2 does not 
directly indicate either hadronic or partonic content of 
the matter. To distinguish between the hadronic and 
partonic content we need more delicate properties of the 
elliptic flow. The calculated results as a rule overesti- 
mate the data. This makes room for the viscosity and/or 
afterburner corrections to the calculated results. In the 
3FD model neither viscosity nor kinetic afterburner at 
the final stage of the expansion were incorporated. 

Elliptic flows of protons and antiprotons at mid- 
rapidity coincide only if these are formed by the single 
fireball fluid, i.e. if baryon-rich fluids (the projectile and 
target ones) are well separated from the fireball fluid in 
rapidity. As seen this is almost the case at 62.4 GeV. 
With the incident energy decrease this separation in ra- 
pidity space becomes smaller, and as a result the the 
proton and antiproton V2 differ from each other more and 
more. This is a consequence of difference in proton and 
antiproton content and in the V2 patters of the baryon- 
rich and fireball fluids. Of course, exceptions from this 
rule are possible if V2 patters of different fluids turn out 
to be very similar by some reason. For instance, this is 
the case for the energy of 19.6 GeV within deconfinement 
scenarios. Typically, V2 of protons and antiprotons be- 
come closer to each other with the incident energy rise, 
similarly to that observed in the experiment. This takes 
place for all considered scenarios because the baryon-rich 
and fireball fluids get more and more separated in rapid- 
ity. 

As for more delicate properties of the elliptic flow men- 
tioned above, contrary to the hadronic scenario, all de- 
confinement scenarios result in proton elliptic flow ex- 
ceeding or being approximately equal (with accuracy of 
the calculation) the antiproton one (with the exception of 
the lowest considered incident energy of 7.7 GeV). The 
latter is in qualitative agreement with the data. The 
lowest energy of 7.7 GeV is already quite close to onset 



of deconfinement within 2-phase-EoS and crossover-EoS 
scenarios [3, [8|. Therefore, the model predicts a kind 
of irregular behavior in this energy region that does not 
agree with the data at present. Within deconfinement 
scenarios the antiproton V2 data are, as a rule, better re- 
produced in the low-momentum region, that is the prime 
region of the applicability of the fluid description. This 
indicates the preference of the deconfinement scenarios. 



III. CONCLUSIONS 

The elliptic flow of produced particles is one of the 
most sensitive observables that brings information about 
the degree of collectivity during the expansion stage of 
heavy-ion collisions. When the collectivity is strong, like 
in the case of ideal hydrodynamics, the elliptic flow takes 
the highest value (the so called hydrodynamic limit). If 
the collectivity is weak, like in a dilute system of weakly 
interacting particles, it is close to zero. The value of 
V2 does not directly indicate either hadronic or partonic 
content of the matter. This is confirmed by approxi- 
mately the same extent of reproduction of the proton 
and antiproton elliptic flow achieved within all scenarios, 
including the purely hadronic one. 

To distinguish between the hadronic and partonic con- 
tent we need more delicate properties of the elliptic flow. 
An indication in favor of the deconfinement scenarios 
is better reproduction of the antiproton elliptic flow in 
the low-momentum region within these scenarios because 
this region is the main domain of the applicability of the 
fluid description. 

In the present version of the 3FD model neither vis- 
cosity nor kinetic afterburner at the final stage of the ex- 
pansion are incorporated. Therefore, it is not surprising 
that the calculated elliptic flow, as a rule, overestimates 
the data. This makes room for the viscosity and/or af- 
terburner corrections to the calculated results. 

The calculated elliptic flows of protons and antipro- 
tons become closer to each other with the incident en- 
ergy rise, similarly to that observed in the experiment. 
In fact, the elliptic flows of protons and antiprotons at 
mid-rapidity coincide only if these are formed by a sin- 
gle fluid, expanding in the center of the colliding system, 
in other words, if this "center" fluid is well separated 
in rapidity from projectile- and target-like leading parti- 
cles, dominantly occupying peripheral rapidity regions. If 
the nuclear stopping rises, as it takes place with incident 
energy decrease, the rapidity gap between the "center" 
fireball and the leading-particles matter shrinks. Then 
the leading-particles matter starts to contribute to the 
mid-rapidity quantities. In its turn, this results in a dif- 
ference between elliptic flows of protons and antiprotons 
that is a consequence of difference in proton and antipro- 
ton content and in the V2 patters of the "center" fireball 
and the leading-particles matter. In the 3FD model this 
mechanism is realized in terms of three interacting fluids: 
the fireball one, consisting of particles newly produced 



near the spacial center of the cohiding system, and two 
baryon-rich fluid {p and t) associated with leading par- 
ticles which traversed the whole system and are finally 
located in longitudinally marginal regions. 

Acknowledgements 

I am grateful to A.S. Khvorostukhin, V.V. Skokov, and 



V.D. Toneev for providing me with the tabulated 2-phase 
and crossover EoS's. The calculations were performed at 
the computer cluster of GSI (Darmstadt). This work was 
supported by The Foundation for Internet Development 
(Moscow) and also partially supported by the Russian 
Ministry of Science and Education grant NS-215.2012.2. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 

[7] 
[8] 

[9] 
[10] 

[11] 
[12] 



S. Voloshin and Y. Zhang, Z. Phys. C 70, 665 (1996) 

hep-ph/94Q7 282 . 

L. Adamczyk et al. [STAR Collaboration], 

larXiv:1301.2348> 

M. M. Aggar wal et al [STAR Collaboration], 

arXiv: 1007.2613 [nucl-ex]. 

Yu. B. Ivanov, V. N. Russ kikh, and V.D. T oneev, Phys. 

Rev. C 73, 044904 (2006) |nucl-th/0503088| . 

V. M. Galitsky and I. N. Mishustin, Sov. J. Nucl. Phys. 

29, 181 (1979). 

A. S. Khvorostukhin, V. V. Skokov, K. Redlich, and V. D. 

Toneev, Eur. P hys. J. C48 , 531 (2006) nucl-th/0605069 . 

Yu. B. Ivanov. larXiv: 1302. 5 766 [nucl-th], to be published 

in Phys. Rev. C. 

Yu. B. Ivanov , Phys. Lett. B721, 123 (2013) 

[arXiv:1211.2579l [hep-ph]]. 

Yu. B. Ivanov. [arXiv: 1304. 1638 1 [nucl-th]. 

V. N. Russkikh and Yu. B. Ivanov, Phys. Rev. C 74 

(2006) 034904 nuc l-th/060 6007 . 

Yu. B. Ivanov, I. N. Mishustin, V. N. Russkikh, 

and L. M. Satar ov, Phys. Rev. C 80, 064904 (2009) 

|arXiv:0907.4"l40l [nucl-th]] . 

Y. .B. Ivanov, Yad. Fiz. 46, 100 (1987) [Sov. J. 

Phys. 46, 63 (1987)]. 



[13] Yu.B. Ivanov, Nucl. Phys. A474, 669 (1987). 

[14] I.N. Mishustin, V.N. Russkikh, and L.M. Satarov, Yad. 

Fiz. 54, 429 (1991) [Sov. J. Nucl. Phys. 54, 260 (1991)]. 
[15] T. Song, S. Plu mari, V. Greco, C. M. Ko and F. Li, 

larXiv:1211.55TT [nucl-th]. 
[16] J. Xu, L. -W. Chen, C. M. Ko and Z. -W . Lin, Phys. 

Rev. C 85, 041901 (2012) |arXiv:1201.3"39n [nucl-th]]. 
[17] J. Steinheimer, V. Koch and M. Bleicher, Phys. Rev. C 

86, 044903 (2012) ,arXiv: 1207.279 1 [nucl-th]]. 
[18] J. C. Dunlop, M. A. Lisa and P. S orensen, Phys. Rev. C 

84, 044914 (2011) ;arXiv: 1107.3"07 8 [hep-ph]]. 
[19] V. Greco, M. Mitrovski a nd G. Torrieri, Phys. Rev. C 

86, 044905 (2012) arXiv: 1201. 4800 [nucl-th]]. 
[20] V. P. Konchakovski, E. L. Bratkovskaya, W. Gassing, 

V. D. Toneev, S. A. Voloshin and V. Voronyuk, Phys. 

Rev. C 85, 044922 (2012) arXiv: 1201. 3320 [nucl-th]]. 
[21] S. A. Bass, M. Belkacem, M. Bleicher, M. Brandstetter, 

L. Bravina, C. Ernst, L. Gerland and M. Hofmann et aL, 

Prog. Part. Nucl. Phys. 41, 255 (1998) nucl-t h/9803035] . 
[22] J. Y. Ollitrault, Phys. Rev. D 46, 229 (1992). 
[23] S. A. Voloshin, A. M. Poskanzer and R. Snellings, 

larXiv:0 809.2949 [nucl-ex].^ 

Nucl. [24] P. Jacobs and G. Cooper, [arXiv:nucl-ex/0008015[ 



